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Narrow line laser cooling is advancing the frontier for experiments ranging from studies of fun- 
damental atomic physics to high precision optical frequency standards. In this paper, we present 
an extensive description of the systems and techniques necessary to realize 689 nm ^5*0 - ^Pi nar- 
row line cooling of atomic ^^Sr. Narrow line cooling and trapping dynamics are also studied in 
detail. By controlling the relative size of the power broadened transition linewidth and the single- 
photon recoil frequency shift, we show that it is possible to continuously bridge the gap between 
semiclassical and quantum mechanical cooling. Novel semiclassical cooling process, some of which 
are intimately linked to gravity, are also explored. Moreover, for laser frequencies tuned above the 
atomic resonance, we demonstrate momentum-space crystals containing up to 26 well defined lattice 
points. Gravitationally assisted cooling is also achieved with blue-detuned light. Theoretically, we 
find the blue detuned dynamics are universal to Doppler limited systems. This paper offers the most 
comprehensive study of narrow line laser cooling to date. 



I. INTRODUCTION 

Narrow line magneto-optical traps (MOTs) are rapidly 
becoming powerful tools in a diverse array of experimen- 
tal studies. These unique and versatile systems have, 
for example, already been used as integral components 
in obtaining fully maximized MOT phase-space densi- 
ties 0, nuclear-spin based sub-Doppler cooling 0, Q, 
all-optical quantum degenerate gases 4], and recoil- free 
spectroscopy of both dipole allowed 5] and doubly- 
forbidden p optical transitions. In the future, narrow 
line MOTs promise to revolutionize the next generation 
of high precision optical frequency standards 7, 8, 9]. 
Narrow line cooling, via the relative size of the transi- 
tion natural width F and the single-photon recoil fre- 
quency shift ur, also displays a unique set of thermal- 
mechanical laser cooling dynamics. In a previous paper 
we explored these behaviors by cooling ^^Sr on the 
^Sq - ^Pi intercombination transition. In the present 
work, we significantly expand upon this discussion and 
provide a comprehensive description of the experimental 
techniques used to realize and study ^Sq - ^Pi cooling 
and trapping. 

Fully understanding narrow line laser cooling requires 
first clarifying the difference between broad and narrow 
Doppler cooling lines. Broad lines, historically used for 
nearly every laser cooling experiment 11], are defined by 
T/lur >> 1. The 461 nm ^^Sr - ^Pi transition shown 
in Fig. 1(a), which typifies a broad line, has, for exam- 
ple, T/ujR ~ 3x10^. In this case, F, or more generally 
the power-broadened linewidth F^;, is the natural energy 
scale. Here, F^; = + s is defined by the saturation 
parameter s = I /Is where / {Is) is the single-beam peak 
intensity (transition saturation intensity). Semiclassical 
physics thus governs the cooling process and the photon 
recoil, although essential to quantitatively understanding 
energy dissipation , serves more as a useful conceptual 
tool than a dominant player in system dynamics. More- 
over, gravity is essentially negligible since the ratio of the 
maximum radiative force to the gravitational force, R = 



hkr/2mg, is typically on the order of 10^, where 27rfi, /c, 
m, and g are Plank's constant, the light field wavevec- 
tor, the atomic mass, and the gravitational acceleration, 
respectively. 

In contrast, narrow Doppler cooling lines are character- 
ized by T/ujR - 1. The 689 nm ^^Sr ^Sq-^Pi transition 
used in this work, for example, has T/cor = 1.6 where 
F/27r {uR/27r = hk^/Aiim) is 7.5 kHz (4.7 kHz). In this 
case, the relevant thermal-mechanical energy scale and 
thus the underlying semiclassical or quantum mechani- 
cal nature of the cooling depends on s. Details of a given 
cooling process are then set by the laser detuning 27:5 = 
A = ujl — uja where ul {^a) is the laser (atomic reso- 
nance) frequency. In particular, S < ^Sq - ^Pi MOT 
dynamics can be divided into three qualitatively distinct 
regimes, hereafter labeled (I - HI), defined by the relative 
size of |A|, F ~ cji^, and F^;. In regime (HI), correspond- 
ing to trapping beam intensities on the order of Is = 
3/iW/cm^ or 5 ~ 1, ur is the natural energy scale. Here, 
single photon recoils and consequently, quantum physics, 
govern trap dynamics. This situation enables limiting 
temperatures of roughly half the photon recoil limit {Tr 
= 2huR/kB = 460 nK, where ks is Boltzmann's con- 
stant) despite the incoherent excitation provided by the 
trapping beams 0, . 

Conversely, for 5 ^ 1 the system evolves toward semi- 
classical physics where Te ^r and hence F^; is the 
dominant energy scale. In this case, cooling and motional 
dynamics are determined by the relative size of F^; and 
A. In regime (II), the A < F^; radiative force produces 
damped harmonic motion and, in analogy to standard 
Doppler cooling, MOT thermodynamics set entirely by 
the velocity dependence of the force 0, . For these 
conditions, the expected S- and 5-dependent tempera- 
ture minima are observed, although with values globally 
smaller than standard Doppler theory predictions. Al- 
ternatively, in regime (I) where A > F^;, the atom-light 
interaction is dominated by single-beam photon scatter- 
ing and trap thermodynamics become intimately linked 
to both the velocity and the spatial dependence of the 
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force. Here, gravity plays an essential role as the ra- 
tio R for the - ^Pi transition is only ~ 16. Con- 
sequently, the atoms sag to vertical positions where the 
Zeeman shift balances 5^ leading to ^-independent equi- 
librium temperatures. 

Narrow line cooling also displays a unique set of S > 
thermal and mechanical dynamics. For these experi- 
ments, the atomic gas is first cooled to /iK temperatures 
and then 5 is suddenly switched from S < to 5 > 0. Sub- 
sequently, the sample evolves from a thermal distribution 
to a discrete set of momentum-space packets whose align- 
ment matches lattice points on a three-dimensional (3D) 
face-centered-cubic crystal 16j. Up to 26 independent 
packets are created with S- and 5-dependent lattice point 
filling factors. Note this surprising behavior occurs in 
the setting of incoherent excitation of a non-degenerate 
thermal cloud. To obtain qualitative insight into the ba- 
sic physics, we begin with an analytic solution to the 
one-dimensional (ID) semiclassical radiative force equa- 
tion. Here, we show that ^ > excitation enables "pos- 
itive feedback" acceleration that efficiently bunches the 
atoms into discrete sets of S- and 5-dependent velocity 
space groups. A simple generalization of the ID model 
is then used to motivate the experimentally observed 3D 
lattice structure. This intuitive picture is then confirmed 
with numerical calculations of the final atomic velocity 
and spatial distributions. Using the numerical calcula- 
tions, we also show that ^ > momentum-space crystals 
are a universal feature of standard Doppler cooling and 
that observations should be possible, although increas- 
ingly impractical, with broad line molasses. 

Finally, we demonstrate that R directly infiuences S 
> thermodynamics, enabling cooling around a S- and 
5-dependent velocity vq where gravity balances the radia- 
tive force. Observed values for vq agree well with numeri- 
cal predictions while cooling is evident in distinctly asym- 
metric cloud spatial distributions that appear in both 
numerical calculations of the cooling process and the ex- 
periment. As with momentum crystal formation, gravi- 
tationally assisted S > cooling is universal to Doppler 
limited systems. In the more typical case where R ~ 10^, 
however, equilibrium temperatures {vq) are on the order 
of hundreds of milli-Kelvin (100 m/s) rather than the 
more useful micro-Kelvin (~ 10 cm/s) values achieved 
with narrow lines. 

The remainder of this paper is organized as follows. 
Section H gives an overview of the 461 nm ^Sq - ^Pi 
MOT used to pre-cool ^^Sr to ^ 2.5 mK. ^P2 metastable 
excited-state magnetic traps that are continuously loaded 
by the - ^Pi cooling cycle are also described. Section 
HI details the highly stabilized 689 nm light source and 
the process used to transfer ^^Sr from the ^5*0 - ^Pi MOT 
to the - ^Pi MOT. Descriptions of the techniques 
used to control and study - ^Pi laser cooling are 
also provided. Sections IV and V then focus on (5 < 
mechanical and thermal dynamics, respectively. Finally, 
section VI explores 5 > cooling and momentum-space 
crystals. Conclusions are given in section VII. 
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FIG. 1: (a) Partial ^^Sr Energy level diagram. Numbers in 
parentheses give Einstein A coefficients (in s~^). Wavelengths 
are in vacuum, (b) Top (x-y plane) view of the Sr cooling 
and trapping apparatus. Blue (red) arrows represent 461 nm 
(689 nm, 679 nm, and 707 nm) trapping (trapping and re- 
pumping) beams. M, mirror; A/4, dual wavelength quarter- 
wave plate; A/2, half- wave plate; DBS, dichroic beamsplitter; 
PBS, polarization beamsplitter; BS, beamsplitter; SMPM, 
single- mode polarization maintaining; ECDL, external- cavity 
diode laser; BA, Bayard- Alpert vacuum gauge; NI, nude vac- 
uum gauge, TSP, Titanium sublimation pump; CC, compen- 
sation coil. 



II. ^So - ^Pi MOT PRE-COOLING 

Pre-cooling ^^Sr to milli-Kelvin temperatures is an es- 
sential requirement for observing - ^Pi cooling and 
trapping dynamics 0. For this purpose, as shown by 
Fig. 1(b), we use a standard six-beam 461 nm ^5*0 - ^Pi 
MOT that is loaded by a Zeeman slowed and transversely 
cooled atomic beam. The atomic beam is generated by 
an effusion oven (2 mm nozzle diameter) whose output is 
angularly filtered by a 3.6 mm diameter aperture located 
19.4 cm from the oven nozzle. Separate heaters maintain 
the oven body (nozzle) at 525 ^C (725 ^C), resulting in a 
measured flux (divergence half-angle) of 3x10^-^ atoms/s 
(19 mrad). The atomic beam is then transversely cooled 
by 2-dimensional 461 nm optical molasses. The elliptical 
cross-section molasses laser beams have a 1/e^ diameter 
of ~ 3 cm (~ 4 mm) along (normal to) the atomic beam 
propagation axis, contain 10-20 mW of power, and are 
detuned from the - ^Pi resonance by -15 MHz. Stray 



magnetic fields in the transverse cooling region are less 
than 1 G. Subsequently, the atomic beam passes through 
a 6.4 mm diameter electro- mechanical shutter and a gate 
valve that allows the oven to be isolated from the rest of 
the vacuum system. 

After exiting the transverse cooling region, the atomic 
beam enters a water cooled 20 cm long constant deceler- 
ation (j~ Zeeman slower 18j with a peak magnetic field 
of ~ 600 G, corresponding to a capture velocity of ~ 
500 m/s. The 461 nm Zeeman slower cooling laser is 
detuned from the - ^Pi resonance by -1030 MHz, 
contains 60 mW of power, and is focused to approxi- 
mately match the atomic beam divergence. The window 
opposite the atomic beam is a z-cut Sapphire optical flat 
that is vacuum sealed via the Kasevich technique 17(], 
broadband anti-reflection coated on the side opposite the 
chamber, and heated to 200 to prevent the formation 
of Sr coatings. The alternating current window heater, 
which produces a small stray magnetic field, is only op- 
erated during the ^Sq - MOT pre-cooling phase. A 
separate compensation coil reduces the slower magnetic 
field magnitude (gradient) to < 100 mG (< 8 mG/cm) at 
the trapping region, located 15 cm from the slower exit. 
Stray magnetic fields at the trap are further nulled by 
three sets of orthogonally oriented Helmholtz pairs. 

The trapping chamber is a cylindrical octagon with 
six 2 3/4" (two 6") ports in the horizontal x-y plane 
(along gravity z-axis). All windows are broadband anti- 
reflection coated. The MOT anti-Helmholtz coils, ori- 
ented such that the axial magnetic field gradient dB^ 
lies along gravity, are mounted on a computer controlled 
precision linear track that allows the coil center to be 
translated from the trapping chamber to a 1 cm x 1 
cm X 4 cm rectangular glass cell. The coils are con- 
structed to provide window-limited optical access to the 
geometric center of the trapping chamber and produce 
axial gradients of 0.819 G/(cm-A). The gradient is linear 
over a spatial range of ~ ± 3 cm in both the axial and 
transverse directions. Current in the coils is regulated by 
a computer controlled servo and monitored with a Hall 
probe. For the ^S^ - MOT, the axial magnetic field 
gradient dB^/dz = dB^ = 50 G/cm. 

The trapping chamber is evacuated by a 40 //s Ion 
pump and a Titanium sublimation pump while the oven 
chamber uses a 40 //s Ion pump. The two chambers are 
separated by a 6.4 mm x 45 mm cylindrical differential 
pumping tube located between the electro-mechanical 
shutter and the gate valve. Typical vacuum levels in the 
oven, trapping, and glass cell chambers during operation 
of the atomic beam are 2xl0~^ Torr, 1.5x10"^ Torr, and 
3xl0~^^ Torr, respectively. 

461 nm cooling and trapping light is produced by fre- 
quency doubling the output from a Ti: Sapphire laser in 
two external buildup cavities that together produce 
> 220 mW of single- mode light. The 461 nm light is then 
offset locked to the - ^P\ resonance by saturated ab- 
sorption feedback to the Ti: Sapphire laser. Relative fre- 
quencies of the transverse cooling, Zeeman slower, and 




Hold Time (sec) Fill Time (sec) 

FIG. 2: Magnetic trap population versus (a) hold time 

for a fixed 2 s loading time and (b) loading time at a hold 
time of 400 ms. Solid lines are exponential fits, td (tl) is the 
measured exponential decay (loading) time. 



trapping laser beams are controlled with acousto-optic 
modulators (AOMs) which are also used as shutters. Ad- 
ditional extinction of 461 nm light is provided by electro- 
mechanical shutters. The intensity stabilized trapping 
beams have 1/e^ diameters of ~ 3 cm, are detuned from 
the ^5*0 - ^P\ resonance by -40 MHz, and typically have 
a total power of 30 mW. For these settings, we find the 
^5*0 - ^P\ MOT population is maximized. Further in- 
creases in, for example, the trapping beam power simply 
increases the cloud temperature. ^5*0 - "^Pi MOTs are 
monitored with a charge-coupled-device (CCD) camera 
and a calibrated photodiode. Typical trap lifetimes, pop- 
ulations, and temperatures are 20 ms, 3x10^, and ~ 2.5 
mK, respectively. 

As shown by Fig. 1(a), operation of the ^S^-^Px MOT 
efficiently populates the ground-state-like met ast able 
excited-state (~ 500 s radiative lifetime 20]) via ^P\ 
^P2 radiative decay. Consequently, - 
MOT lifetimes are typically limited to 10 - 50 ms 
|2l|, [23, 23] . To overcome this loss process, ^P2 popula- 
tion is re-pumped to the ground-state via the ^P2 
^Si ^Pi ^So channel by driving the 707 nm ^P2 - 
^Si and 679 nm ^Pq - ^Si transitions with two external 
cavity diode lasers (ECDLs). Each laser is locked to a ref- 
erence cavity that is simultaneously locked to a frequency 
stabilized helium neon laser. Double-passed AOMs are 
then used to tune the absolute laser frequencies. After 
passing through a single-mode polarization-maintaining 
optical fiber, the co-propagating 707 nm and 679 nm laser 
beams are expanded to a 1/e^ diameter of ~ 1 cm and 
delivered to the ^Sq - ^Pi MOT. An AOM located before 
the beam expansion optics allows for rapidly turning the 
beams either on or off (< 1 /is transition time). At the 
trap, the 707 nm (679 nm) beam contains 1.5 mW (2.5 
mW) of power, resulting in an optical re-pumping time 
of < 100 /is. With both lasers operating, the - ^Pi 
MOT population and lifetime are typically enhanced by 
10 X and 15 x, respectively, with the former value limited 
by atomic beam induced trap loss. 

Along with contributing to ^5'o - ^Pi MOT loss, ^Pi 
^P2 radiative decay continuously loads milli- 
Kelvin atoms in the ^P2(m=l,2) states into a magnetic 
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FIG. 3: 689 nm laser system. Solid (dashed) lines indicate op- 
tical beams (electrical connections). A/4, quarter- wave plate; 
A/2, half- wave plate; M, mirror; CM, curved mirror; PBS, 
polarization beamsplitter; BS, beamsplitter; 01, optical iso- 
lator; SMPM, single-mode polarization maintaining; ECDL, 
external- cavity diode laser; DL, diode laser; TSVI, tem- 
perature stabilized and vibration isolated; PD, photodiode; 
AOM, acousto-optic modulator; EOM, electro-optic modula- 
tor; DBM, double balanced mixer; INT, integrator; AM, am- 
plitude modulation input; CUR (PZT), current (piezo-electric 
mounted grating) modulation input. 



trap formed by the - MOT quadrupole mag- 
netic field 111 111 il lil HI Importantly, these 
samples are expected to display a wealth of binary col- 
lision resonances that arise due to an interplay between 
anisotropic quadrupole interactions and the local mag- 
netic field m 1^1 l3lj . Qualitatively similar processes 
and hence, collision resonances, are predicted for polar 
molecules immersed in electrostatic fields 32] . Studies of 
metastable Sr collision dynamics will thus likely impact 
the understanding of a diverse range of physical systems. 
To pursue these studies, we plan to first load state 
atoms into the qua drup ole magnetic trap and then me- 
chanically translate 1^3 ^he sample to the glass cell cham- 
ber. Subsequently, a tight loffe-Pritchard magnetic trap 
will be used to perform a variety of collision experiments. 

As a first step in this direction we have loaded > 10^ 
^P2(m=l,2) state atoms into the quadrupole magnetic 
trap and achieved trap lifetimes > 1 s. For these mea- 
surements, the trap is first loaded by operating the 
^S'o - ^Pi MOT for a variable fill time. The 461 nm cool- 
ing and trapping lasers and the atomic beam are then 
switched off. Following a variable hold time, the 461 nm 
trapping beams and the 707 nm and 679 nm re-pumping 
beams are switched on, enabling the magnetic trap 
population Nm to be determined from 461 nm fluores- 
cence 25j. Figure 2(a) shows Nm versus hold time for 
a fill time of 2 s while Fig. 2(b) gives Nm versus fill 
time at a fixed hold time of 400 ms. For both, dB^ = 50 
G/cm, giving a magnetic trap depth of ~ 30 mK. The 
- ^Pi MOT lifetime (steady-state population) is 16 
ms (2.8x10''). The observed 1.23(4) s exponential decay 
time agrees well with the expected ~ 1.3 s vacuum limited 
value extracted from measurements reported in Ref. 27] 
at higher pressures. In contrast, the significantly shorter 



0.55(2) s fill time implies that additional loss processes 
are operative during the loading phase, possibly due to 
interactions between atoms in the or states with 
^5*0 - ^Pi MOT atoms or the atomic beam 34]. This idea 
is supported by the magnetic trap loading rate. For the 
parameters used here, the observed rate of 2.7(5) x 10^ 
atoms/s is ~ 2x smaller than the theoretically predicted 
value of 5.4(1.8) X 10^ atoms/s 24] where the theory only 
accounts for relevant branching ratios in the ^Pi 

^P2(m=l,2) radiative cascade. Similar discrepancies 
between observed and predicted loading rates were re- 
ported in Ref. 



III. ^5^0 - ^Pi MOT LOADING AND DETECTION 

Exploring 689 nm - ^Pi narrow line cooling dynam- 
ics requires a laser system whose short-term linewidth is 
small compared to the 7.5 kHz transition natural width. 
In addition, the absolute laser frequency must be refer- 
enced, with similar stability, to the ^Sq - ^Pi transition 
and be tunable over a range of ~ 10 MHz. Figure 3 shows 
the 689 nm laser stabilization and control system consist- 
ing of a master-slave ECDL, a temperature stabilized and 
vibration isolated passive optical reference cavity, and a 
Sr saturated absorption spectrometer. 

The linewidth of the master ECDL is first narrowed by 
locking the laser to a stable optical reference cavity via 
the Pound-Drever-Hall technique 1^3 • The cavity con- 
sists of high reflectivity Zerodour substrate mirrors that 
are optically contacted to a Zerodour spacer. The mea- 
sured cavity finesse (free spectral range) at 689 nm is 
- 3800 (488.9 MHz), giving a linewidth for the TEMqo 
mode of ~ 130 kHz. To isolate the cavity from environ- 
mental perturbations, the cavity is suspended by two thin 
wires inside a temperature stabilized can that is evacu- 
ated to < 10~^ Torr and mounted on vibration damp- 
ing material. The absolute cavity frequency is tuned 
by double-passing the 689 nm light sent to the cavity 
through a 1 GHz AOM. The electronic feedback, with an 
overall bandwidth of ~ 2 MHz is divided into a slow loop 
that adjusts the piezo-electric mounted ECDL grating 
and a fast loop that couples to the diode laser current. 
With the cavity lock engaged, in-loop analysis shows that 
jitter in the cavity-laser lock is < 1 Hz. 

To further evaluate the performance of the 689 nm 
laser system, the short-term laser linewidth is deter- 
mined by beating the cavity-locked 689 nm light against 
a femto-second comb that is locked to a second optical 
cavity 36]. Here, we find a short-term linewidth of < 300 
Hz, where the measurement is limited by the optical fiber 
connecting the 689 nm light to the femto-second comb. 
Next, the absolute stability of the 689 nm laser is eval- 
uated by beating the cavity locked 689 nm light against 
a femto-second comb that is locked to a Hydrogen maser 
via a fiber optical link to NIST |37l. From these measure- 
ments, the laser-cavity system drifts ~ 400 mHz/s and 
has a 1 s stability of < 4x10"^^ (i.e., < 180 Hz) with the 
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FIG. 4: ^5'o - ^Pi MOT timing diagram. Images at top 
show the atomic cloud at various stages of the coohng and 
compression process. From left to right, the first four images 
show in-situ images of the cloud while the last frame shows 
the cloud after 25 ms of free expansion from the S = -520 
kHz, s = 75 single- frequency MOT. N, cloud population; Tm, 
cloud temperature. 



former value limited by cavity drift and the latter value 
limited by the effective maser noise floor. To eliminate 
the slow 689 nm cavity induced drift, the master ECDL 
is next locked to the - ^Pi resonance via saturated 
absorption feedback to the 1 GHz AOM. For added sta- 
bility, a DC magnetic field is applied to the Sr vapor cell 
and the spectrometer is set to perform frequency mod- 
ulation spectroscopy on the ^6*0 - ^Pi(m=0) transition. 
With the system fully locked, the 1 s stability (drift rate) 
is then - 4x10"^^ (< 80 mHz/s). 

A portion of the 689 nm master ECDL output is next 
used to injection lock a 689 nm slave diode laser. The 
slave laser output, after double-passing through an AOM 
used for frequency shifting and intensity chopping, is then 
coupled into a single-mode polarization-maintaining op- 
tical fiber. Upon exiting the fiber, the 689 nm light, 
containing up to 6 mW of power, is expanded to a 1/e^ 
diameter of 5.2 mm and divided into three equal intensity 
trapping beams. Dichroic beamsplitters are then used to 
co-align the 689 nm and 461 nm trapping beams. The 
trapping beam waveplates are 3 A/4 at 461 nm and A/4 
at 689 nm. 

- ^Pi cooling and trapping dynamics are monitored 
either by in-situ or time-of-flight (TOF) fluorescence im- 
ages collected with an intensified CCD camera. The cam- 
era is set to view the cloud in either the horizontal x- 
y plane or nearly along gravity. The x-y plane (along 
gravity) images have a spatial resolution of 21 /im/pixel 
(37 /im/pixel). For in-situ images, the 461 nm trapping 



beams are pulsed on for 10 - 50 /is immediately after 
the atoms are released from the trap while for TOF im- 
ages, the atoms are allowed to first freely expand for a 
variable amount of time. We have verified that in-situ 
images recorded with 461 nm pulses are identical to di- 
rect images of the in-trap 689 nm fluorescence aside from 
an improved signal-to- noise ratio. Typical TOF flight 
times are 20 - 35 ms. To determine cloud temperatures, 
gaussian fits are performed to both the in-situ and TOF 
images. The temperature Tm is then given by Tm = 
(m/4/cBtF^)(KF - r^)^ where tp is the flight time and 
Rp (r) is the TOF (in-situ) 1/e^ radius of the cloud. 

Figure 4 depicts the ^5*0 - ^Pi MOT loading procedure 
Q. As outlined above, ~ 3x10^ atoms are first pre- 
cooled to - 2.5 mK in a ^Sq - ^Pi MOT. (Note that 
for the remainder of this paper, the - ^Pi MOT 
is loaded without the 707 nm and 679 nm re-pumping 
lasers.) At time t = 0, the 461 nm light and the atomic 
beam shutter are switched off, dB^ is rapidly lowered to 
3 G/cm, and red-detuned, broadband frequency modu- 
lated 689 nm trapping beams are turned on. 10 ms later 
and for the following 50 ms, the cloud is compressed by 
linearly increasing dB^ to 10 G/cm. Frequency modula- 
tion parameters for the 689 nm trapping beams are set 
to provide complete spectral coverage of the ^5*0 - ^Pi 
MOT Doppler profile and, as shown below, manipulate 
the cloud size at the end of the magnetic field ramp. 
Subsequently, at t = 60 ms, the frequency modulation is 
turned off and the atoms are held in a single-frequency 
MOT. Overall, as shown by the images at the top of Fig. 
4, this process reduces the sample temperature by more 
than three orders of magnitude while only reducing the 
cloud population by a typical factor of 3 - 4, giving final 
temperatures of ~ 1 /iK and populations of ~ 10^. Typ- 
ical single frequency trap lifetimes and spatial densities 
are ~ 1 s and ~ 5x10^^ cm~^, respectively. 



IV. ^ < MECHANICAL DYNAMICS 

As outlined in section I, ^ < ^6'o - ^Pi MOTs provide 
a unique opportunity to explore three qualitatively dis- 
tinct laser cooling regimes whose underlying mechanics 
are governed by either semiclassical or quantum mechan- 
ical physics. In this and the following section, we will 
describe the unique experimental signatures for regimes 
(I) - (HI) and give detailed explanations for the observed 
trapped atom behavior. 

Insight into regime (I) and (II) thermal-mechanical dy- 
namics is provided by the semiclassical radiative force 
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FIG. 5: (a) Semiclassical radiative force versus position (bot- 
tom axis, Vx = Vy = 0) and velocity (upper axis, x = y = 0) 
for a range of detunings. Corresponding (b) trap potential 
energy in the z direction and (c) in-situ images of the ^5*0 - 
^Pi MOT. Dashed lines in (c) are calculated maximum force 
contours. For each, s = 248. 



FIG. 6: (a) Semiclassical radiative force versus position (bot- 
tom axis, Vx = Vy = 0) and velocity (upper axis, x = y = 0) 
for a range of intensities. Corresponding (b) trap potential 
energy in the z direction and (c) in-situ images of the ^So - 
^Pi MOT. Dashed lines in (c) are calculated maximum force 
contours. For each, S = -520 kHz. 



where x = {x^y^z}, dB = {dBx, dBy^ dBz}, and fi = 
{qjIIb/K) where gj = 1.5 (/i^) is the ^Pi state Lande 
g-factor (Bohr magneton), s' > s accounts, along a sin- 
gle axis, for saturation induced by the remaining four 
trapping beams. Figure 5(a) shows the -^So - ^Pi radia- 
tive force dit s = s' = 248 for a range of 5 values while 
Fig. 6(a) shows the force at (5 = -520 kHz for a range 
of s = s' values. As described in section I, the quali- 
tative nature of the force and hence the resulting trap 
mechanical dynamics depends on the relative size of A 
and F^;. For regime (I), corresponding to \S\ > 120 kHz 
in Fig. 5(a) or s < 248 in Fig. 6(a), A > F^; and the 
3D radiative force acts only along a thin shell volume 
marking the outer trap boundary. Here, the trap bound- 
ary roughly corresponds to positions where the radiative 
force is peaked. This situation, as shown by Figs. 5(b) 
and 6(b), produces a box potential with a gravitation- 
ally induced z-axis tilt. Hence, in the x-y plane, motion 
consists of free-flight between hard wall boundaries while 
along the z-axis, mechanical dynamics are set by the rel- 
ative size of the radiative force "kicks," gravity, and the 
cloud thermal energy. As shown in section V, the thermal 
energy is small compared to the gravitational potential 
energy. Moreover, the ratio of the maximum radiative 
force to the gravitational force, R ^ 16. Thus, the atoms 
sink to the bottom of the trap where they interact, along 
the z-axis, with only the upward propagating trapping 
beam. 

As S decreases in Fig. 5(a) or s increases in Fig. 6(a), 
the trap mechanically evolves to regime (II) where A < 
F^; produces a linear restoring force and hence, damped 
harmonic motion ITEj . Consequently the trap po- 
tential energy assumes the U-shaped form familiar from 
standard broad line Doppler cooling. As the trap moves 



more fully into regime (II), perturbations to the poten- 
tial energy due to gravity become less pronounced. One 
expects, therefore, that the cloud aspect ratio will evolve 
toward the 2:1 value set by the quadrupole magnetic field. 

The intuitive descriptions developed above are directly 
confirmed by Figs. 5(c) and 6(c) which show in-situ im- 
ages of the ^5*0 - ^Pi MOT along with overlaid maxi- 
mum force contours calculated from Eq. ^ For excita- 
tion conditions corresponding to regime (II), the cloud 
approaches the 2:1 quadrupole magnetic field aspect ra- 
tio. In contrast, for regime (I) the cloud x-y width is 
determined largely by the separation between x-y force 
maxima or alternatively, by the wall separation for the 
x-y potential energy box. In the vertical direction, the 
atoms sink to the bottom of the trap where the lower 
cloud boundary zq is defined by the location of the z-axis 
potential energy minima which is, in turn, proportional 
to the position where the Zeeman shift matches the laser 
detuning. As S increases, zq shifts vertically downward, 
an effect predicted in Fig. 5(b) and clearly revealed in 
Fig. 5(c). To quantify this relationship. Fig. 7 shows 
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FIG. 7: Vertical cloud position versus ^ for s = 248 and dB^ 
= 10 G/cm. The solid line is a linear fit. 
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FIG. 8: Modematching broadband cooled atoms to the singl( 
frequency MOT. (a) Broadband modulation spectra. Dashe 
vertical lines give the effective detuning corresponding to tt 
measured cloud size, (b) In-situ images of the broadban 
cooled atoms. For each, the temperature is 8.5(0.8) /iP 
Dashed lines give the maximum force contour for the S ■ 
-520 kHz, s = 75 single-frequency MOT shown in the bo 
tom frame, (c) Broadband to single-frequency MOT transf( 
efficiency versus broadband spectra number. 



zq versus S along with a linear fit to the data givin 
dzo/d(^ = 27T/{jiidBz) = 0.509(4) /im/kHz, in agreemer 
at the 5% level with the expected linear slope of 0.478(/ 
/im/kHz. 

Unique - ^Pi MOT mechanical dynamics are also 
manifest in the transfer between the broadband and 
single-frequency cooling stages. Here, it is important to 
note that following broadband cooling, the cloud verti- 
cal and horizontal 1/e^ radii (r^ and r^, respectively) are 
set by the final value for dB;^ and the spectral separa- 
tion between the blue edge of the modulation spectrum 
and the - ^Pi resonance. Figure 8(a) shows three 
typical modulation spectra while Fig. 8(b) shows in-situ 
images of the corresponding broadband cooled clouds at 
the end of the magnetic field ramp. For each, the mea- 
sured temperature is 8.5(0.8) jaK. Overlaid dashed lines 
give the maximum force contour for the subsequent 6 = 
-520 kHz, 5 = 75 single-frequency MOT shown in the 
bottom frame. Here, dynamics similar to Fig. 5 are 
observed: as the broadband modulation moves closer to 
resonance, the cloud density distribution becomes more 
symmetric and the aspect ratio evolves toward 2:1. Fig- 
ure 8(c) shows the broadband to single- frequency MOT 
transfer efficiency TE versus broadband spectra number. 
Clearly, TE increases as the overlap between the broad- 
band cooled cloud and the single-frequency MOT force 
contour increases, indicating that TE can be optimized 
by "mode-matching" r^ and r^ to the single-frequency 
MOT box. This idea is supported by open circles in Fig. 
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dB^ = 10 G/cm. 



8(c) which give predicted TE values based on the relative 
volume of the broadband cooled cloud and the single- 
frequency radiative force curve. Overall, the prediction 
reproduces observed TE values with a discrepancy for 
spectra 1 and 2 that likely arises from the broadband 
MOT density distribution. 

Further evidence for broadband to single-frequency 
mode-matching is provided by the single-frequency MOT 
population N versus S shown in Fig. 9. Here, solid cir- 
cles give N versus S when r^ = 0.79 mm and r^ = 1.56 
mm. According to the Fig. 5 model, the single-frequency 
MOT acquires these dimensions at 5 = -1637 kHz, which 
lies within 2% of the S value where N is peaked. Hence, 
we again find that optimal transfer occurs when r^ and r^ 
are matched to the single-frequency radiative force con- 
tour. Additional confirmation for this effect is provided 
by open circles in the figure which give N when r^ and 
Th are optimized at each 6. To understand the trend in 
this case toward smaller N as 5 decreases, note the non- 
zero slope for the blue edge of the modulation spectrum 
sets the minimum effective detuning Sm that can be used 
during the broadband cooling phase. In analogy to Fig. 
5, Sm sets minimum values for iz and r/j,. Thus mode- 
matching becomes progressively more difficult as 5 de- 
creases and the single-frequency radiative force contour 
shrinks, leading to the observed decrease in N. 



V. ^ < THERMODYNAMICS 

S < ^ So - ^Pi MOTs display a rich variety of ther- 
modynamic behaviors that are directly linked to the me- 
chanical dynamics explored in section IV. Figure 10(a) 
shows the MOT equilibrium temperature Tm versus S for 
saturation parameters ranging from s = 7.5 to s = 1800. 
For large S and 5, corresponding to regime (I), A > F^; 

r, Tm is basically ^-independent. Insight into this be- 
havior is provided by Fig. 11(a), which details the unique 
regime (I) connection between trap thermodynamics, the 
spatial dependence of the radiative force, and the relative 
size of the radiative force and gravity. Recall that for A 
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FIG. 10: ^So - ^Pi MOT temperature versus (a) detuning 6 
for various intensities and (b) intensity s at ^ = -520 kHz. In 
(a), DT is standard Doppler theory. In (b), solid circles are 
experimental data while the solid and dotted lines give the 
Doppler limit (hVE/'^kB) and the single-photon recoil limit 
{2hujR/kB), respectively, (c) Global scaling factor applied 
to Doppler theory in order to match regime (II) data versus 
intensity. 



> r^;, the cloud sags to the bottom of the trap where in- 
teractions occur, along the vertical z-axis, with only the 
upward propagating trapping beam. Moreover, due to 
polarization considerations and the free-flight motion ex- 
ecuted by atoms in the horizontal plane, horizontal beam 
absorption rates are more than 4x smaller than the verti- 
cal rate. Trap thermodynamics, therefore, are dominated 
by a balance between gravity and the radiative force due 
to the upward propagating beam. Here, it is important 
to realize that as S changes, the z-axis atomic position zq 
self adjusts such that the effective detuning, A - /idB^^o, 
remains constant. Consequently the trap damping and 
diffusion coefficients, and thus the equilibrium tempera- 
ture, remain constant. 

To obtain a quantitative expression for Tm under these 
conditions, we first find the damping coefficient a by Tay- 
lor expanding 



hkv 

mg. 



4(A - kvz - iidB^zf 



Y2 



(2) 



where {dF/dv) = a is evaluated dit Vz = 0^ z = zq. 
Solving F(0, zo) = and using R = {hkV /2mg), the 5- 



FIG. 11: (a) Schematic depiction of regime (I) cooling. 
Atoms sag to the bottom of the potential energy (PE) curve 
where interactions occur primarily with the upward traveling 
trapping beam. The inset shows the excitation geometry, (b) 
Numerical factor Ni? in Eq. |S] versus intensity. 



independent effective detuning is given by 



A - jidBzZ \/Rs - s' -1 



(3) 



which, in combination with Eq. [SJ states that the scat- 
tering rate depends only on i?, 5, and T. Substituting 
this expression into {dF/dv)^ we obtain the damping co- 
efficient 



2hk^VRs -s' -1 



(4) 



Next, the diffusion coefficient Dp is calculated by substi- 
tuting Eq. Olinto the single-beam scattering rate. Dp is 
then given by 



D^ 



2R 



(5) 



Combining Eqs. 01 and O the predicted (^-independent 
equilibrium temperature is 



R 



a| 2^R- s'/s-l/s 



2kB 



Nr 



(6) 



where we have used ^/s ~ y/s -\- 1 for 5 ^ 1. As shown 
by Fig. 11(b), the numerical factor Nr is approximately 
2 over the entire relevant experimental range. Regime 
(I) temperatures, therefore, depend only on F^;. To test 
this prediction. Fig. 10(b) shows Tm versus s for a fixed 
large detuning 6 = -520 kHz. For the central portion of 
the plot where regime (I) dynamics are relevant, we find 
good agreement with Eq. ^ 

For small 6 in Fig. 10(a), the system evolves from 
regime (I) to regime (II) where A < Te^ ^ ^- As this 
transition occurs, trap dynamics change from free-flight 
to damped harmonic motion. Here, one expects thermo- 
dynamics similar to ordinary Doppler cooling including 
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FIG. 12: The approach to thermal equihbrium at ^ = -520 
kHz. (a) Measured temperature versus time for s = 248. The 
inset gives a semi- log plot for times greater than 5 ms. The 
solid line is a linear fit. T, instantaneous temperature; To, fi- 
nal equilibrated temperature (b) Measured exponential equili- 
bration time versus s. The dashed, dotted, and solid lines are 
predictions from standard Doppler theory, the single-beam 
theory outUned in the text, and quantum theory, respectively. 



5- and 5-dependent minima with equilibrium values hav- 
ing the functional form l^, 



T{s,A) 
To 



To(4|A|/rs)- 

{hrE)/{2kB) 



-1 (l + 4(A/r£;) 



(7) 



where To, realized at A = r£;/2, is a generalized version 
of the 5 <C 1 Doppler limit. As shown by the solid lines 
in Fig. 10(a), Eq. Q correctly reproduces the functional 
shape of the data. As shown by Fig. 10(c), however, 
matching the absolute data values requires multiplying 
Eq. [3 by a s-dependent global scaling factor (< 1) whose 
value decreases with 5, leading to temperatures well be- 
low the standard Doppler limit Tq. In contrast to ordi- 
nary Doppler cooling, the cloud thermal energy in regime 

(II) is thus not limited by half the effective energy width 
of the cooling transition. Notably, we find this surprising 
result cannot be explained either by analytic treatments 
of Eq. Qor semiclassical Monte-Carlo simulations of the 
cooling process. The Monte-Carlo simulations, in fact, 
simply reproduce standard Doppler theory. Finally, as s 
approaches unity in Fig. 10(b), the trap enters regime 

(III) where F ~ cor ~ ksT/h and thus the cooling be- 
comes fully quantum mechanical. Here, we obtain a min- 
imum temperature of 250(20) nK, in good agreement 
with the quantum mechanically predicted value of half 
the photon recoil temperature Tji/2 = hcjR/kB = 230 
nK 14]. 

The approach to thermal equilibrium also displays sig- 
natures of regimes (I) - (III). Figure 12(a) gives the cloud 
temperature versus time for a (5 = -520 kHz, s = 248 
single-frequency MOT. As shown by the figure inset, the 
data is well fit by a single exponential aside from the 
rapid decrease for times < 5 ms that arises from atoms 
falling outside the trap capture velocity |3, |^. Recall 
from Fig. 10 that scanning s at this detuning provides 
access to all three cooling regimes, each of which is char- 



acterized by a unique evolution toward thermal equilib- 
rium. Specifically, for regimes (I) and (II), semiclassi- 
cal Doppler theory predicts and equilibration time r = 
m/2\a\ where a is given in regime (I) by Eq. E| and, 
according to ordinary Doppler theory, in regime (II) by 



Ahe lAi 



a = 



1 + s + 4(A/r)2 



(8) 



Conversely, quantum theory predicts that for regime 
(III), T follows Ji,] 



(9) 



where is the initial cloud temperature. As a test of 
these predictions. Fig. 12(b) shows r versus s at (5 = - 
520 kHz along with r values predicted by the above three 
theories. Note that for each s value we find the tem- 
perature versus time is well fit by a single exponential 
aside from the ~ 5 ms long rapid decrease typified by 
Fig. 12(a). In regime (II) at s = 4520, r = 14(2) ms, in 
good agreement with Doppler theory which gives r ~ 12 
ms. As expected for regime (I), corresponding to 10 < 
s < 4000, r follows equilibration times predicted by the 
single-beam damping coefficient while in regime (HI) at s 
= 6.1, r = 3.1(4) ms, consistent with the Eq. [Ul quantum 
mechanical predictions. 



VI. ^ > COOLING AND MOMENTUM-SPACE 
CRYSTALS 

Tuning to ^ > during the single-frequency cooling 
stage reveals two fundamental and unique physical pro- 
cesses: (1) the creation of well-defined momentum pack- 
ets whose velocity space alignment mimics lattice points 
on a face-centered cubic crystal and (2) laser cooling 
around a velocity vq where the radiative force balances 
gravity. In the following, we explore these two effects in 
detail. 

Basic insight into ^ > momentum packet formation 
can be obtained by considering the elementary problem 
of ID atomic motion in the presence of two counter- 
propagating S > light fields. For simplicity, we assume 
dB = 0. According to Eq. ^ an atom with initial veloc- 
ity Vi will preferentially interact with the field for which 
^ • ^ > 0. Hence, the absorption process preferentially 
accelerates rather than decelerates the atom, further de- 
creasing the probability for absorption events that slow 
the atomic motion and enabling "positive feedback" in 
velocity space. This process terminates for final veloci- 
ties Vf satisfying 



t;/ » i ( A ■ 



(10) 



which has a linear A dependence and scales with F^;. 
Neighbor atoms with initial velocities around vl undergo 
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FIG. 13: One dimensional velocity bunching due to single- 
beam ^ > absorption for (a) ^ = 80 kHz, and (b) 5 — 200 
kHz. Solid, short dash, and dashed lines correspond to initial 
velocities of 0.1 cm/s, 5 cm/s, and 10 cm/s, respectively, (c) 
Three-dimensional ^ > momentum space structure. Final 
momenta corresponding to the black, gray, and white circles 
arise from three-, two-, and single-beam interactions, respec- 
tively. 



similar "positive feedback" acceleration and ultimately 
achieve final velocities near Vf. To quantify this latter 
effect, we first simplify the problem by neglecting the 
beam for which k ■ Vi < 0. The equation of motion for 
the atomic velocity is then 



dv 

m 



hkV 
~2rn 



1 



4(A - kvY 

f2 



(11) 



which can be solved analytically for the interaction time 
t as a function of the velocity v 



t{v,Vi) 



- 4Akv' + ((1 + s)T' + 4A')v 

-k^f - 4Akv^ + ((1 + s)T^ + 4.A^)vi] (12) 
o 



Figures 13(a) and 13(b) show plots of Eq. [Hat J = 80 
kHz and 6 = 200 kHz, respectively. For both, 5 = 30 or 
s = 300. The chosen Vi values of 0.1 cm/s, 5 cm/s, and 
10 cm/s span the velocity distribution for a 12 /iK cloud. 
Assuming an interaction time of 25 ms. Fig. 13(a) clearly 
shows that due to laser induced acceleration, the entire 
range of initial velocities is rapidly bunched to a signif- 
icantly reduced range of final velocities. Considering a 
fully ID situation, this result implies the cloud is divided 
into two well defined and oppositely moving packets. As 
5 increases to 200 kHz in Fig. 13(b), the 5 = 30 atom- 
light interaction becomes sufficiently weak for ~ that 
two final velocity groups centered around Vf ^ S cm/s 
and Vf ^ 2b cm/s appear. In contrast, we find that for s 
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FIG. 14: Top view in-situ images for a range of S and s at 
\di3\ — and t// = t\^ = 25 ms. Arrows in the uppermost 
center frame show the propagation direction for the horizontal 
molasses beams. The initial cloud temperature is 11.5(2) /xK. 
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FIG. 16: Broad line momentum space crystals. Calculated 
horizontal (a) final versus initial velocity and (b) correspond- 
ing spatial distributions for 461 nm ^ > optical molasses. 
For each, in = 500 /is and the initial temperature is 2.5 mK. 



FIG. 15: Calculated horizontal ^ > final versus initial ve- 
locity for (a) various ^ at s = 30 and (b) various s at ^ = 120 
kHz. (c) and (d) give the corresponding spatial distributions. 
For each, t// = 25 ms and the in = temperature is 11.5 /xK. 



= 300 the transition to two groups does not occur until 
S ~ 400 kHz. Accounting again for the full ID symme- 
try, the 5 = 30 case corresponds to dividing the t = 
cloud into three groups, two that move in opposite direc- 
tions with relatively large final velocities and one that, 
in comparison, is nearly stationary. Overall, Eq. ITTlthus 
predicts that for J > the atomic cloud evolves into a 
discrete set of momentum-space packets. As s decreases 
at fixed the number of packets increases while, as pre- 
dicted by Eq. El the mean velocity for a given packet 
scales with both S and s. 

Generalizing this analysis to the full 3D molasses beam 
geometry leads to the structure shown in Fig. 13(c): a 
3D array of momentum-space groups whose alignment 
mimics the lattice points on a face-centered-cubic crys- 
tal. From symmetry considerations, cube corners corre- 
spond to three beam interactions in which the cloud is 
divided into two pieces along each coordinate axis. From 
Figs. 13(a) and 13(b), these points have a 6- and 5- 
dependent mean velocity and appear at relatively small 
detunings. For larger 5 values, two-beam interactions fill 
points along the Fig. 13(c) corner connecting lines. In 
analogy to 5 = 30 in Fig. 13(b), atoms in these points 
have Vi ^ along a single axis and thus remain nearly sta- 
tionary. Along the remaining two axes, however, > 0, 
enabling acceleration to larger Vf. The 3-beam to 2-beam 
transition, as shown in Figs. 13(a) and 13(b), occurs at 
progressively larger 6 values as s increases. Together, the 
two processes form a total of 20 divided groups with 8, 4, 
and 8 packets in the top, middle, and bottom layers of the 
cube, respectively. Finally, as 6 increases further, atoms 
with Vi ^ along two axes experience acceleration only 
along a single axis, producing the I-beam lattice points 
shown as 6 open circles in the Fig. 13(c) cube face cen- 
ters. 



Figure 14 shows an array of top view (slightly off ver- 
tical) in-situ cloud images for intensities ranging from s 
= 30 to 5 = 1040 and detunings spanning = 80 kHz 
to S = 200 kHz. For each, \dB\ = and the atom-light 
interaction time is fixed at in = 25 ms (ty = 25 ms) in 
the horizontal x-y plane (along z-axis) molasses beams. 
The initial in = W = cloud temperature is 11.5(2) 
fiK. Here, the observed cloud evolution agrees well with 
the qualitative predictions developed from the Fig. 13 
model. As 5 increases at fixed s, sets of n-beam lattice 
points sequentially fill with n = 3 filling first followed by 
n = 2,1. Moreover, as s increases, transitions between 
the n-beam processes occur at progressively larger 6 val- 
ues. Finally, for fixed (5, the mean lattice point velocity, 
proportional to the resulting lattice point spacing in the 
Fig. 14 images, scales with s as expected from both Eq. 
Eland Fig. 13. Note that qualitatively similar dynam- 
ics are observed for \dB\ ^ 0. The detailed changes to 
the cloud evolution that result from this situation are 
described below. 

The intuitive understanding of ^ > dynamics devel- 
oped above is confirmed by Fig. 15. Here, we show nu- 
merical calculations based on Eq. Qof the final horizontal 
velocity and spatial distributions for the 5 = 30 column 
and (5 = 120 kHz row in Fig. 14. For the calculations, 
\dB\ = 0, ti:/ = 25 ms, and the initial cloud temperature 
is 11.5 /iK. The Figs. 15(c) and 15(d) spatial distribu- 
tions should be compared to cube lines in the x-y plane 
along the Fig. 14 x-y molasses beam propagation direc- 
tions. Importantly, this fully ID model reproduces both 
the Fig. 14 observations and the Fig. 13 predictions for 
the ^- and 5-dependent lattice point filling factors and 
mean lattice point velocity. As expected, the temper- 
ature of each packet in its moving frame is lower than 
the tj:/ = ty = ms atomic cloud, a result arising from 
the previously described velocity bunching effects and di- 
rectly connected to cloud shape asymmetries observed in 
both the experiment (note the sharp outer cloud edges) 
and the Fig. 15 theory. Notably, however, only two ver- 
tical layers are observed in Fig. 14 while the Fig. 13 
model predicts three. As explained below, this appar- 
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FIG. 17: Side view in-situ images for s = 30, \dB\ = ( 
= 140 kHz, Ih = 25 ms, and a range of iv times. The 
= tv = cloud temperature is 11.5(2) /xK. Calculated 
final versus initial velocity for the same set of parameters. 
Corresponding calculated spatial distributions. 



ent discrepancy arises from novel gravitationally indu 
z-axis dynamics. 

At this point, it is important to realize that Eq. □ 
is semiclassical and thus does not account for dynam- 
ics that depend on the relative size of T and lor. The 
5 > momentum crystals observed here, therefore, are 
a universal feature of Doppler limited systems. Similar 
processes should then occur with broad lines such as the 
461 nm - ^Pi transition. To test this possibility. Fig. 
16 shows numerically simulated final horizontal velocity 
and spatial distributions for a 2.5 mK ^^Sr cloud excited 
by 461 nm ^ > optical molasses. For the calculation, 
\dB\ =0,tH = 500 /is, and 5 = 2. As clearly indicated by 
the figure, structures very similar to those shown in Fig. 
14 can be generated over length scales consistent with 461 
nm cooling beam diameters and typical - ^Pi MOT 
temperatures. Moreover, we find that by increasing ei- 
ther 5 or ti^, structures with contrast identical to those 
shown in Figs. 14 and 15 can be created. Note the model 
does not account for spontaneous emission induced ran- 
dom walk heating, a process the would tend to smear the 
contrast between individual momentum packets. This 
omission, however, should not significantly affect Fig. 16 
since random walk heating scales as VT 12] while, from 
Eq. n the directional acceleration scales as F. Thus, as F 
increases from the narrow to the broad line case, random 
walk heating becomes progressively less important. 

Although the physics underlying S > momentum- 
space crystals is thus fully operative for broad lines, de- 
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FIG. 18: Measured vertical position of the upward moving 
layer versus ty at 6 = 140 kHz and s — 30. The solid line is 
a fit to the model described in the text. 



tails of the crystal formation process will make exper- 
imental observations difficult. From Eq. El ^/ foi" a 
given packet is proportional to F. Hence, broad line 
packet velocities are orders of magnitude larger than 
those achieved in the narrow line case. In Fig. 16, for 
example, Vf ~ ± 60 m/s. Visualizing the entire structure 
then requires imaging light with an optical bandwidth of 
~ 260 MHz or an equivalently broad optical transition. 
This situation should be compared to the experiments 
performed here where Vf < 1 m/s, and hence sufficient 
imaging bandwidth is obtained with the F/27r = 32 MHz 
^5*0 - ^Pi transition. In addition, F sets the lattice point 
spacing, or alternatively the molasses beam diameter, re- 
quired to achieve a given contrast between individual lat- 
tice points. For the contrast shown in Figs. 14 and 15, 
for example, the required broad line molasses beam diam- 
eters grow to tens of centimeters, making experimental 
observations impractical. 

Figure 17(a) shows in-situ images of the 689 nm S > 
molasses when the cloud is viewed in the horizontal x-y 
plane at 45^ to the x-y axes. For the images, \dB\ = 0, 5 
= 30, 5 = 140 kHz, in = 25 ms, and ty is varied from ty 
= ms to ty = 25 ms. Here, the apparent contradiction 
between Fig. 14 and the Fig. 13(c) prediction for the 
number of vertical layers is shown to occur due to gravity 
induced dynamics. Increasing ty from ty = ms to 
ty = 6 ms, for example, creates three vertical layers, 
as predicted by Fig. 13(c). For ty > 6 ms, however, 
the lower two layers slowly merge together, becoming a 
single cloud along the vertical direction for ty > 21 ms. 
To understand this process, recall that the central layer 
in the Fig. 13(c) cube corresponds to atoms with near 
zero z-axis velocities. For the chosen S and in the absence 
of gravity, therefore, these atoms would remain near Vz = 
0. As shown by the numerical calculations in Figs. 17(b) 
and 17(c), however, gravity accelerates the central layer 
into resonance with the downward propagating molasses 
beam, causing the two downward propagating layers to 
merge. For the ty = 25 ms time used in Fig. 14, this 
process is complete. Hence, only two layers are observed 
with the more (less) intense packets in the Fig. 14 images 
corresponding to the lower two (uppermost) cube layers. 

Comparing the numerically simulated S > spatial dis- 



tributions with Figs. 14 and 17(a) reveals that theoreti- 
cahy predicted \dB\ = lattice spacings are ~ 2x larger 
than observed. This result occurs due to stray magnetic 
fields. Here, an independently measured ~ 100 mG/cm 
permanent chamber magnetization spatially shifts the ef- 
fective detuning as the atoms move outward, causing an 
apparent deceleration and thus reduced lattice spacing. 
This effect can be seen most clearly by measuring the ver- 
tical position of the upward moving cube layer zt versus 
ty. For the upward moving layer, Eq. ^predicts that the 
atoms are accelerated to a velocity vq where the radiative 
force balances gravity. The cloud then moves upward at 
vq and hence, experiences an effective acceleration = 
0. To test for a magnetic field induced non-zero ^e, Fig. 
18 shows Zt versus ty for s = 30 and S = 140 kHz. The 
solid line in the figure is a fit to the simple kinematic 
equation zt = zq -\- voty - (^e^y)/2. Here, vq is treated 
as an initial velocity since for s = 30 and 5 = 140 kHz, 
the atoms are accelerated to 'Uq in < 2 ms. From the 
fit, we find = -0.98(16) m/s^, consistent with a stray 
gradient dB^ = 100 mG/cm. Once this gradient is in- 
cluded in numerical calculations of Eq. the measured 
and calculated lattice point spacings agree. 

Finally, the exceptionally sharp velocity and thus spa- 
tial distributions for the upward propagating layer in Fig. 
17 imply velocity compression beyond the bunching ef- 
fects described earlier. In fact, as shown by Fig. 19(a) 
which depicts the composite gravitational and radiative 
force for the upward moving atoms, stable cooling occurs 
around vq where the composite force is zero. From Eq. 
12 '^0 is given by 

^o = l[^+^Vs{R-l)-l^ (13) 

which depends linearly on A and scales approximately as 
r^;. Fig. 19(b) shows measured values for vq versus S at 
5 = 10 and s = 75, clearly demonstrating the expected 
linear ^-dependence. To obtain quantitative comparisons 
with theory, we next perform linear fits to the data and 
then obtain predicted values for vq from numerical cal- 
culations that include the 100 mG/cm field gradient dis- 
cussed above. From the experiment, we find {dvo/dS) = 
644(17) /im/(kHz-s) at s = 10 and 634(19) /im/(kHz-s) 
at s = 75. These slopes agree at the 10% level with the 
predicted values of {dvo/dd) = 572(5) /im/(kHz-s) at s = 
10 and 638(20) /im/(kHz-s) at s = 75. Moreover, exper- 
imentally observed absolute values for vq are reproduced 
by the calculations at the level of 20%, in good agreement 
with the expected Eq. intensity dependence. 

The equilibrium temperature of the upward moving 
layer can be calculated by following the same procedure 
used to derive Eq. [SI Here, however, we set jidBzZ = 
in Eq. [3 and perform the Taylor expansion about Vz = 
Vq. The (5-independent effective detuning is then 

A-kvz ^ VRs -s' -1 

r ~ 2 ^ ^ 




(cm/s) Detuning (kHz) 



FIG. 19: ^ > cooling, (a) Composite gravitational and S 
> radiative force versus Vz- Arrows show the direction for 
force induced velocity changes. Stable cooling occurs around 
a velocity vq, labeled by a circle, where gravity cancels the 
radiative force, (b) Measured vo versus 5 for s = 10 and s = 
75. Solid lines are linear fits. 



which leads to damping and diffusion coefficients iden- 
tical to Eq. El and O respectively. Hence, the expected 
equilibrium temperature is given by Eq. Thus, grav- 
ity, via the ratio again plays an important role in 
narrow line thermodynamics, in this case enabling S > 
cooling. Unfortunately, this prediction cannot be ac- 
curately verified due to spatial overlap among the hor- 
izontal plane packets in the upward moving cube layer. 
The observed and theoretically predicted sharp vertical 
spatial distributions in Figs. 17(a) and 17(c), however, 
strongly suggest that S > cooling is operative in the 
experiment. Finally, note that although this same cool- 
ing mechanism should occur for broad lines where R ~ 
10^, ^0 in this case has impractical values on the order of 
100 m/s. Moreover equilibrium temperatures are large, 
at roughly W0x{hTE/2kB) 200 mK for the - ^Pi 
transition even at s = 1. 



VII. CONCLUSIONS 

In summary, narrow line laser cooling exhibits a wealth 
of behaviors ranging from novel semiclassical dynamics 
wherein gravity can play an essential role to quantum me- 
chanically dominated sub-photon recoil cooling. In the 
^ < semiclassical case, trap dynamics are set by either 
hard wall boundaries or a linear restoring force. Qualita- 
tive differences between these two situations are reflected 
in both the atomic motion and the equilibrium thermo- 
dynamics. Here, mechanical dynamics range from free- 
flight in a box potential to damped harmonic oscillation. 
Accordingly, equilibrium temperatures range from detun- 
ing independent values scaled by the power-broadened 
transition linewidth to detuning dependent minima well 
below the standard Doppler limit. As the saturation pa- 
rameter approaches unity, the trap enters a quantum me- 
chanical regime where temperatures fall below the pho- 
ton recoil limit despite the incoherent trapping beam ex- 
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citation. For ^ > 0, the cloud divides into momentum- 
space crystals containing up to 26 well defined lattice 
points and the system exhibits ^ > gravitationally as- 
sisted cooling. These surprising S > behaviors, which 
again occur due to an incoherent process, are theoreti- 
cally universal features of Doppler limited systems. Ob- 
servations should therefore be possible, although difficult, 
with broad line optical molasses. Perhaps a similar S > 



mechanical evolution also occurs for atoms, such as the 
more typically employed Alkali metals, that support both 
Doppler and sub-Doppler cooling. 
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